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ABSTRACT 

Among the tools available for the study of the dark energy driving the expansion of the Universe, Baryon Acoustic Oscillations (BAO) 
and their effects on the matter power spectrum are particularly attractive. It was recently proposed to study these oscillations by map- 
ping the 21cm emission of the neutral hydrogen in the redshift range 0.5 < z < 3. We discuss here the precision of such measurements 
using radio-interferometers consisting of arrays of dishes or north-south oriented cylinders. We then discuss the resulting uncertainties 
on the BAO scales and the sensitivity to the parameters of the Dark Energy equation of state. 
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c3 1. Introduction 



' One of the goals of current cosmological research is to fully 

characterize the "dark energy" that drives the apparent accelera- 

i tion of the expansion of the Universe. Among the available tools 

■ JAlbrecht et al.ll2006l) . much interest has been generated by the 
' features in the matter power spectrum that result from Baryon 

Acoustic Oscillations (BAO) in the pre-recombination Universe. 

These oscillations cause a peak in the matter correlation function 

, at a comoving distance, s ~ 1 05 /i^'Mpc, equal to the acoustic 

■ horizon at recombination. Equivalently, the BAO induce "wig- 
gles" in the matter power spectrum with peaks at comoving wave 
numbers k - {n + 1/2)7t/s for n = 2,4, The theoretical po- 
sitions of these peaks are known to a precision of order 1 % so 
they can be used as reliable standard rulers to study the expan- 
sion history. 

After first being seen in the C MB anisotropy spectrum 
' dMauskopf et"ani2000l: iHinshaw et al...2008.) . BAO effects were 
subsequently observed in the low-redshift (z < 1) galaxy cor- 
relation function and power spectrum by the SDSS survey 
(Eisens tein et al. 2005; Percival et al. 2007) and the 2dGFRS 
survey dCole et al]l2005l) . With future data in the redshift range 
0.5 < z < 3, it is hoped to provide precise constraints on dark- 
energy parameters. It has been propos ed to do this with r edshift 
surveys using galaxy optical spectra (Bas set et al J 120051) or HI 
emission ( Abdala & Raw l ings 2005) or by mapping in three di- 
mensions the pat tern of Lyman-a absorption of distant quasars 
(IWolfetalJl2005h . 

An alterna ti ve ele gant approach was recently proposed by 
iPeterson et alj (l2006h and further developed by IChang et alj 
(|2008|) . They propose to use neutral hydrogen (HI) as a tracer 
of matter and simply map out the 21cm emission with an an- 
gular resolution that is insufficient to detect individual galax- 
ies (where most HI is concentrated) but sufficient to observe the 
BAO wiggles. At redshift z = 1.5, the acoustic horizon subtends 
0.0334rad on the sky which can be resolved with telescopes of 



size of order 100m. These would be considerably smaller and 
cheaper than the 1km elements necessary to de tect individual 
high redshift galaxies (Abdala & Rawlings ;2005h . 

Mapping the matter distribution using HI 21 cm emis- 
sion as a tracer has been e x tensively di scussed in literature 
(e.g. iFurlanetto eFal] (l2006h ; iTegmark & Zaldarriaga (2008)). 
Several proje cts, such as LOFAR (Rottgerin g et alj 1200 6) or 
MWA (.Bowman et alj 12007.) aim at detecting the reionization 
epic (z ~ 10). Detecting BAO features a round z ~ 1 using H I 
radio emission has also been discussed bvl Wvithe et al.l (l2007h . 

In Section |2] we review the expectations for the HI power 
spectrum PniiK) and its relation to the large-scale structure mat- 
ter power spectrum PissiK). Section [3] describes the interfero- 
metric observations that are assumed for this study and Section 
ID presents a simplified procedure to use these observations to 
reconstruct the HI power spectrum in a small volume near the 
zenith. While this procedure is not entirely realistic, it allows 
us to give estimates of the reconstructed noise power spectrum, 
Pnoiik) (Fig.|5]). Section|5]discusses the problem of subtraction of 
foreground and background radio sources. Finally, in Section |6] 
we describe the constraints on the cosmological parameters that 
can be derived from the determination of the BAO peaks in the 
power spectrum. 

2. The HI power spectrum 

The mean HI brightness temperatur e (assumed to be mu ch 
greater than the CMB temperature) is dBarkana & Loebir2007h 



Th,(z) = 0.0466m/:- 



(l+zffHi(z) 



^Hiho 



(1) 



[Q.A + il+z)^Q.M]^'^ 3.5x10-4 

where Q.hi ~ 3.5 x 10^ dZwaan et al.ll2005h is the current mean 
cosmological HI density and 



fmiz) 



(HI/H),=Q 



(2) 
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is the fraction of hydrogen in atomic form relative to this fra ction 
at z = 0. Studies of Lyman-or absorption dWolf et al]|2005 ') indi- 
cate the fniiz = 1.5) ~ 3 implying fmiz = 1.5) = 0.00033/:. 
This is three orders of magnitude below the brightness of ex- 
tragalactic radio source emission, ~ 34007y 5r"', coiTesponding 
to frs ~ 03K. However, these sources, as well as Milky Way 
emission and CMB radiation, have smooth frequency distribu- 
tions, allowing them, in principle, to be subtracted. 

We are interested in the spatial variations of the HI temper- 
ature about the mean. Consider a small cube of volume V at 
mean redshift z coiTesponding to a HI frequency v - vo/(l + z). 
The cube covers a solid angle AQ = Aa x A6 {a = i.a., 6 - 
dec.~ 0) and redshift range Az corresponding to a frequency 
range Av = voAz/(l + zf'- The present volume of the cube is 



V = didnA^ Az 



(c///o)4AQAz 



(3) 



where for a flat ACDM universe, the comoving angular distance 
dr and Hubble distance dn are 



dr - 



{c/Ho)dz 



^JnA + ^Mil+zy 



dn = 



cIHq 



VnA+nM(l+z)3 



Inside the cube, the sky brightness can be expanded using 
functions that satisfy periodic boundary conditions. For small 
cubes (Aa, AS, Az <K 1) the functions can be taken to be com- 
plex exponentials of {a, 6, z): 

THi{a,5,z) = I^yr/,;(n)e2-X«»-/A»+«.5/Afi+H,./Az) ^4^ 

where the Tniin - («„, ns, nS)), are defined for integer «„, rts, «-. 
The wave vector k is related to n by 



K/y 



lima 



h - 



2Km 



djAa djAS 

The HI power spectrum is: 

PHiik) = <|r„,(^)|2) 



Inn, 
1^^ - L 

'' dHAz 



(6) 



the average being over the modes k with wave numbers near k. 
The number of such modes in the cube and in the interval A^ is; 



A^AA- 



Vk^ Ak _ (c/Ho)dT(zfAQ.Azk^ Ak 



: 7.25 X lO"* 



dT(z)f AQAz / 



27r/5 \ 0.075/1 Mpc-i 

1/2 



c/Ho 

xl^A + il+zf^MJ 



Ak/k 
0.2 



(7) 



Since the |r(A:)p are random number with variance equal to P{k) 
the precision with which P(k) can be measured is P{k)/ '\/Nm- 

The HI power spectrum is expected to be similar to the 
galaxy power spectrum. Figure |5] shows the expected galaxy 
spectru m normalized to agre e with that measured by SDSS at 
z ~ 0.1 dTegmark et al.ll2004 J). Also shown is this spectrum ex- 
trapolated to z = 1.5 assuming Q.m,^a = 0.27,0.73. The spec- 
trum has BAO peaks at 

k„^(n + 1/2)7t/s ^ 0.075/iMpc'' " ^ ^^^ n = 2,4,6 (8) 

where s - 105/? 'Mpc is the sound horizon at recombination. 
The relative crest to trough amplitudes at the first three peaks are 
expected to be AP/P ~ 0.13, AP/P ~ 0.09, and AP/P ~ 0.05. 



3. Interferometric Observations 

Observations of GHz radiation can be performed with interfer- 
ometers consisting of arrays of reflectors and receivers. In this 
paper, we consider two types of arrays. The first, drawn schemat- 
ically in Fig. [T] consist of north-south oriented cylindrical re- 
flectors pointing towards the zenith. The reflectors have dipole 
receiving antennae deployed along their axes. The information 
to be extracted from a pair of such cylinders is illustrated in Fig. 
[1] The amplitude for each receiver is sampled over a time f,„, 
after which a Fast-Fourier Transform (FFT) is performed to sep- 
arate frequency components. We take f,„, to be sufficiently long 
(~ 35 fisec) to give excellent resolution on the BAO features in 
the radial direction. The receivers within a given cylinder are 
then combined by a FFT to form beams covering pixels in dec- 
lination of width ~ A/Dg for a wavelength A and cylinder length 
Dg. This requires that receivers be spaced by A/2 in order to 
avoid large side-lobe contamination. The number of receivers 
per cylinder is thus of order 500 for Dg - 100m. The amplitudes 
from two cylinders separated in uA in the east-west direction 
can than be correlated to form a "visibility", V(a(t), 6, v, u) that 
is defined for discrete values of declination, 6, and frequency, v, 
defined by the two FFT's. The right ascension, a{t) is defined by 
the time of observation, t. 



S(a(t),6,v) 




FFT 

N-S 



(5) L 



FFT 

35|J,sec 



FFT 

N-S 



Lt;^)-^ V(a(t), 5,v.u) 



Fig. 1. Two cylinders instrumented with dipole receiving an- 
tennae along their axes. After a sampling time f,„, (taken here 
to be 35fisec), an FFT is performed for each receiver to sepa- 
rate frequency, v, components. Receivers within a given cylinder 
are then combined by (north-south) FFT to form beams in dec- 
lination, 6. The amplitudes from the two cylinders can then be 
combined to form visibilities, V(a{t), 6, v, u) where a{t) is the 
right-ascension determined by the time of the observation and 
where the cylinder separation isuA. 



The second type of array considered in this paper consists 
of dishes the are orientable in declination. The readout for pairs 
of dishes is the same as for the cylinders in Figure [1] except that 
there is no north-south FFT so that the declination 6 is defined 
by the common dish pointing. Additionally, the dish separation 
has components in the north-south and east-west directions so u 
is a vector 

The cylinders of Fig. [1] give high resolution information in 
the declination and frequency directions. In order to provide high 
quality information in the right-ascension direction, it is nec- 
essary to have several cylinders with a range of reflector sepa- 
rations Au. Schematics of the three arrays considered here are 
shown in Figure|2l Two arrays, a and b, are aiTays of fixed cylin- 
drical reflectors of width D„ = 10m and length Dg - 100m. The 
east-west configuration of the five cylinders in configuration a is 
chosen to give uniform sensitivity over the k range necessary for 
BAO studies. As such, they are placed at positions that are inte- 
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a n 



100m 



10m 



b '' 



100m 



c 



10m 



5m 



Fig. 2. The three telescope configurations considered: a: "un- 
packed cylinders"; b: "packed cylinders"; c: "packed dishes". 
In the unpacked array, the five cylinders are separated from 
the most western cylinder by distances of « x 10m with n - 
0,4,5,7,13. 




Xu„sina' 



Fig. 3. The pointing configuration of telescope pairs (at the ter- 
restrial equator for simplicity). The left panel shows two cylin- 
ders separated in the east-west direction and pointing toward the 
zenith. The right panel shows two dishes separated in the north- 
south direction and pointing in a direction 6 from the zenith. 



gral multiples of the cylinder width as shown in the figure. On 
the other hand, the ten cylinder in configuration b are adjacent 
to form a "packed array" giving redundant information over the 
interesting k range. 

The third system, c, is an array of D„ xDg - 5m x 5m dishes 
(taken to be squares for mathematical convenience). Each dish 
is instrumented with one receiver. The field-of-view in both an- 
gular directions is ~ A/ Da- In order to survey the sky, the dishes 
must be pointable in declination. As for the cylinder arrays, the 
right-ascension is determined by the time of observation. 

As illustrated in Fig. [1] each readout-FFT sequence gener- 
ates a set of amplitudes Sadfc indexed by the four integers: a 
giving the right ascension, a„, determined by the time of obser- 
vation; (d, f) denoting the set of amplitudes generated by FFT 
coiTesponding to pixels in declination-frequency space centered 
on 6d, Vf (only one d per readout for dishes); and c defining the 
reflector (dish or cylinder). The amplitudes are the sum of signal 
and noise: 



The angular and frequency response are determined by the 
lobe functions, LgdjcW ,d' ,v'), peaked ?A a' - 6' - v' - and 
normalized so that the expectation value of the squared signal is 

{\S%f) = fda'd6'dv'L(a',S', v')T(a' + a„ 6' + 6d. v' + v/) 

For the rest of this note we make the approximation that Ladfc is 
independent of {a, d, f, c). The noise contribution is normalized 
so that the mean signal-to-noise is given by f/Tsys 



(\S"adff} = T,yAda' dd' dV Lia', 6', v') 



(10) 



The visibility for direction (a, d) and frequency / for cylin- 
ders (or dishes) c and c' separated in space by uAf is: 



Vadfu - SadfcS gclfc' 



(11) 



The noise-noise contribution is a random number of vanishing 
mean and variance determined by r^y^: 



(12) 



{WZ%\^f^ = T^y^^da' d5' dV L{a',6', V) 

Following Fig. |3] the signal-signal contribution to the visibility 
has an expectation value 

<^fld/«> = jda'd6'dv'L(a', 6', v') T(a' + a„, 6' + 5d, V -H v/) 

X exp [-27r/(MQ, sin a' H- m^ sin(5rf H- <5'))] (13) 
For a' ,d' « 1 and 5^/ <K 1 this simplifies to 

I da' dd' dv L{a', 6', v') T(a' + a^, 6' + 6d, V + v/)e"^'"'"'^ 
where (6»' =(a',5'))- 



4. Reconstruction of HI Power spectrum 

The visibilities, Yadfa can be combined to give information on 
the HI distribution. For a given readout, a Fourier transform of 
Yadfu over u yields a map of the field of view. We are more inter- 
ested in the power spectrum which is related to the Fourier trans- 
form of Yadfu over right-ascension, declination and frequency: 

V(k,U) = i y -KdjKaMTh&,^dnk,z,)y IniuA (15) 

NdNfNa j^. ^ 

where there are Nd xN/xNa pixels densely and uniformly cov- 
ering the cube (Aa, A6, Av). In the absence of noise, substituting 
the expansion (|4]i into ( fT4l i shows that V{k, u) is proportional to 
the product of Tniik) and the Fourier transform of the lobe func- 
tion. More generally, the visibilities are due to noise and signal 
so we have 



Vik,u) = 



f r(Jt) Fik, u) 



i da'd5'dv'Ua',S',v') (16) 



" adfc 



c not , c sig 
'^ adfc ^ '^ adfc 



(9) 



where r(k) - rHiik)+r„oi{k) (plus foreground/background con- 
tributions) and where the dimensionless "form factor" for the 
mode k is 

r da' dd' dv'L(a', 6', y')e'i''<'H(i+z)-y'eiffidTk-2nu) 

F(k,u) = ^ . (17) 

j da'dd'dv'L{a',6',v') 
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If L is symmetric, F is a real function. 

The estimate of r(^) based on the visibility of one reflector 
pair with separation u is denoted r(^, «). From ( fT6b . it is given 
by 



r(A:,H) =C 



V{k,u) 
F(k,u) 



(18) 



where the "calibration constant", C, relates the power spectrum 
to the measurements and form factors: 



C = 



Vv 



T jda'd6'dv'L(a',6',v') 



(19) 



Since more than one value of u can be used to estimate r(^), 
one can make a weighted average of the r(^, u) derived with 
( fTSl l. In the next subsection, we will show that the noise power 
is inversely proportional to F(k, u) so it is reasonable to have 
weights that are increasing functions of \F(k, u)\ ~ F(k, u). It is 
most convenient to take the weights to be equal to F{k, u): 



T{k) 



Y.„F(k,u)r{k,u) ^_i:„m,«) 



(20) 



where the sum is over those values of u with F(k, u) i^Q (Other 
choices of weighting may give slightly lower noise power but at 
the price on introducing wiggles in the noise power spectrum.) 

4. 1 . The noise 

The noise contribution to the V(k,u) ( fTsT i is the mean of 
NaNdNf random numbers of variance given by (fTZt . This gives 

(y„oi(k, m)) = and 



1/2 



I^V„oi{k,ut) 



2\I/2 



^jNaNdNf 



(21) 



The noise, r„„,(^, u) reconstructed using (fTsT l is inversely pro- 
portional to F(k, u) with an expectation value, Pq, given by 



Po = {\r„,Ak,u)F(k,ut) =^ 



f2 NdNfNa 



(22) 



The survey volume to pixel ratio V/NaNdNf depends on the 
interferometer configuration. In all cases the number of pixels 
along the frequency direction is Nf - f,„/Av where ?,„, is the in- 
tegrating time between electronic readouts. In the angular direc- 
tions, cylinders and dishes differ. Dishes have only one angular 
beam, so we have NgNd = tio,/ti„, where f,o, is the total obser- 
vation time of the survey. Cylinders, with FFT beam-forming, 
have Na - tt„t/ti„t and A^^ = A6/{A/Ds) where Dg is the cylinder 
length. Combining with (O this gives 



V 



d^duAQ. 



(1 + z)^NaNdNf 



ttofVO 

d^dnAn A/Ds 



(dishes) 



t,o,vo A6 



(cylinders) 



(23) 



The values of Pq thus differ for cylinders and dishes. For 
cylinders we have 



/ Aa/t,o, \ ^JnA + ^Md + z)^ I ^Hiho 



\27t/3 X 10^ sec j (I + z)^ fHi(z)^ \3.5 x IQ-^ 



For dishes we have a much larger Pq because of the smaller sky 
coverage: 



Pq = 4.6 X W(h-'Mpcy 



; / Tsys \ I dr 



\50KJ \c/Hq 

AaiUo, \ V^A + nM(i+z)3/ ^„,h '-^ 



(25) 



^271/3x10^ sec j (I + z)- fHi(z)- \3.5 x IQ-^y 
Using (|20] |. the expectation value of the noise for the mode 



k is 



<ir„.,(^)p) = 



PoNu(k) 



iZu\F(k,u)\f N,(k){\F(k,u)\)l 



(26) 



where Pq is given by ( |24] | or ( |25] ) and Nu(k) is the number of vis- 
ibilities available for estimating r(^). The noise is thus inversely 
proportional to N„(k){F)^. We will see that the large value of f o 
for dishes compared to cylinders will be compensated for by the 
larger value of A^„. 

4.2. The form factor 

To estimate the form factor, we use the lobe function 



\2/ \2/ \2 

smna'Dal A\ (smnd' Dg/A\ /sin7rv'f,„( ] 



L(a',6',v') 



For a cylinder, D„ is the width, Dg is the length, and f,„, is the 
integrating time. This lobe function is a reasonable approxima- 
tion both for the diffraction limited beam in right-ascension and 
for the FFT formed beams in declination and frequency (if high 
frequency components are filtered out). The form factor for this 
lobe is a product of triangle functions A(x) (A(0) = 1, A(x) - 
for \x\ > 1). For ?,„, > 30jU.sec, the frequency lobe is sufficiently 
naiTow that its Fourier transform is near unity for interesting val- 
ues of k. In this case, the form factor is the product of the two 
angular factors: 



F(k,u) = A 



dTka/ln — Ua 
DJA 



A 



drks/ln — Us 

dJa 



(27) 



For cylinder arrays, us - 0. The most important characteristic of 
this form factor is that it vanishes as ^ — > (i.e. k <K InDa/Adj). 
This is because only correlations between separated pairs of 
cylinders are used (no self-correlations) which makes the tech- 
nique insensitive to noise drifts but at the price of insensitiv- 
ity to structure on high angular scales. Because of the trian- 
gular form of the form factor, for the unpacked cylinder array, 
there are generally two values of m^ that have non-vanishing 
F(k, u) for a given k: the ones with Ua just above and just be- 
low ka lyinDa I (Adj)}. For the other two arrays, there are more 
than two combinations. 

Figure |4] plots A^„(^) and 2„ F(k,u) averaged over orienta- 
tion of k for the three telescope configurations. The dish aiTay 
has many more pairs for a given k than the cylinder aiTay and 
this works to compensate for its higher noise. 

4.3. The power spectrum 

The power spectrum P(k) - Piss + Pnoi can be found by calcu- 
lating a weighted average of the |r(^)p for k in an interval Ak: 



(24) P(k) 



Y^k W(k)\T(k)\^ 
Y.kmk) 



(28) 



^ 

n^ 

w = 



O 






O 
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Fig. 4. 2« F(k, u) (solid lines) and Nu(k) (dotted lines) both av- 
eraged over the orientation of k. The top pair of curves is for the 
dish array (configuration c in Figure |2]i, the middle pair for the 
packed cylinder array (configuration b in Figure |2]i, the bottom 
pair for the unpacked cylinder array (configuration a in Figure 
|2]l. The vertical lines give the positions of the first four BAO 
peaks. 

where r(^) is given by ( |20] l. A reasonable choice is W{k) - 
(Yju F(I^' "))^ in which case 



i:k\i:uF{k,u)nk,u)\' 

ZkiZuFik,u))^ 



P(k) 



Using ( |26] | we find a simple expression for the noise: 



(29) 



^noi\'^) 



{{I.uF{k,u))^) 



where the averages are over k. As expected from ( l26t . P„„i(k) is 
inversely proportional to Nu{F}^. 

Figure |5] shows the P„oi(k) calculated in the three configura- 
tions for a four-month observation of 2;7rsr at z = 1.5 (Az - 0.2). 
The noise for the packed aiTay is considerably smaller than that 
for the unpacked array because of the higher value of A^„. The 
larger value of Pq for the dishes is compensated by the larger 
value of A^„ so the noise is comparable to that for the cylinder 
arrays. Note that the dishes accomplish this with fewer electron- 
ics channels, 400, than the 2500 (5000) channels needed for the 
unpacked (packed) cylinders. 

For all configurations, the noise diverges as fc — > because of 
the filtering by the lobe, resulting in decrease in the form factor 
and an increase in P„oi{k). The minimum k with full sensitivity 
is 



2nD,„i,JA 
driz) 



0.038/1 Mpc' 



Dmin/A 3.lGpC 

10m/0.525m dr 



(31) 



where D„„„ is the minimum cylinder or dish separation. For the 
numerical values of dj and A in this formula we have taken 




k(hMpc-0 



Fig. 5. The P(k) and uncertainties for a four-month observation 
of 2nsr at z = 1.5 (Az - 0.2). The solid lines labeled P„oi give 
the noise power calculated with ( l30l l for the three configurations 
of Figure |2] The dashed lines shows the SDSS P{k) at z ~ 
and extrapolated to z = 1.5. The dotted line shows the peak-to- 
peak amplitude of the BAOs at z - 1.5. The solid lines labeled 
CTp shows the uncertainty in the measured power spectrum for 
Ak/k = 0.2 calculated with ( l32b for configuration a and b. The 
vertical lines show the positions of the first four BAO peaks. 

z = 1.5 {Q.M,^A - 0.27,0.73) and the D,„,„ value for config- 
urations b and c of Fig.|2] Figure |6] plots k,„i„(z) for D,„,„ - 10m 
and D,„in - 5m. To have good sensitivity at the first BAO peak, 
10m elements are acceptable at z = 1.5 but smaller elements are 
needed at z = 0.5. 

The uncertainty, crp, in the measured power spectrum aver- 
aged over an interval Ak is 



{wikfink)]") 



1/2 



crp = ^ , , „, (32) 



WTkiWik)) 

where N^k is the number of modes in range Ak given by (|7]i. If 
the noise dominates the power spectrum, this is 

PolNuikfy^ 
(30) crp = „_:,. [ ,.,, (33) 



'Nrk{(ZuF(k,um 

Figure|5]shows the resolution crp calculated for a four-month 
observation of Insr at z = 1 .5 (Az - 0.2) for the two cylinder 
configurations a and b, the dish configuration giving intermedi- 
ate results. The resolution is less than the BAO amplitude shown 
in the figure. For the packed array b, P„oi < Plss at the first BAO 
peak so no improvement can be made by increasing the observ- 
ing time. For the unpacked array c, increasing the observing time 
by a factor 10 would leak to a resolution similar to that for the 
packed array. 

;. The value 



liPnoi > Plss, ctp is of order (Po/Nuik))/ ■ 



of this quantity for unpacked cylinders and calculated at the first 
BAO peak is plotted as a function of redshift in Fig. |6] Also 
shown is the peak-to-peak amplitude of the BAO oscillation. The 
resolution relative to the BAO amplitude rises slowly with red- 
shift. 

5. Radio sources 

Extragalactic radio sources and Milky Way Synchrotron radia- 
tion have power law frequency spectra that allow them to be sub- 



R. Ansari et al.: Reconstruction of HI power spectra with radio-interferometers to study dark energy 




104 



I 



BAO amplitude 



100- 



10- 



P /n"' 

Ak 



1 



Fig. 6. The top panel shows the z dependence of k„,i„ - 
InD/iAdr) for D = 10m (solid line) and D - 5m (dashed line). 
The two horizontal lines show the position of the first two BAO 
peaks. In the bottom panel, the solid line shows Pq/ y/NAk for an 
exposure of 2nyr - sr cylinders in configm^ation a. The dashed 
line shows the BAO peak-to-peak amplitude at the first peak 
(lower panel). 



tracted, in principle, from the HI spectrum. We consider here ra- 
diation for extragalactic sources because their random positions 
on the sky lead to a flux highly dependent on angular position. 
The radio source contribution to the sky brightness (01 is 



Trs(a, 6, v) 



y^i 



2kbV^ ^ 



6{a - ai)5(5 - 5i) 



where the sum is over radio sources / of angular position (a,, 5,), 
\A2GHz flux (in Jy) O,-, and spectral index is y,. The radio 
source power spectrum is determined by the Fourier transform: 



r„(A:) = 



d\dH{\+zf 
Thi Vv vo 



/ 



X I dadSdvTfsict, 6, v)e 



-i^djk^a+drhd+dtiil+zfk.vlva) 



The V integral can be done by developing v '' to first order around 
the mean frequency in the box, v - vo/(l + z). The result is 



r„(A:) = -i(-\)"exp(-id,k-J{\+z)) 



di Av „^, 



k^y/v VQ 



2kbvlTHi 



rjika, ks) 



where n = n^ is defined by (|4]i, y ~ 0.7 is the mean spectral 
index (iBinnev & Merrifield.1998.) , and 77 is a sum over the radio 
sources in the angular region AorAJ: 



77(fe„, h) = (r + 2)-' Yj-ji + 2)(1 + zy--^%e 



-idrikoa+ksS) 



(36) 



(Note: z is the redshift of the box, not of the radio source). The 
power in the mode k is then 



ir„(^)p = 






kldn 

(l+^)2(r+i) 



Az(r + 2)2 



2kbvlTHi) 



AQ-'\7]ika,ket 



(37) 



This spectrum is highly anisotropic because fluxes are correlated 
in frequency but not in angular direction. 
The expectation value of I77P is 

AQ.-' {m,,ks)f) = J ^O-^O^ ~ 67.7 7^,^-1 (33) 

where we use the spectrum from Ijacksoni (1200 4)). Half of the 
integral comes from the ~ 5000 sources per steradian with fluxes 
greater than 0.06Jy. This gives 

"^ - \l/2 



[aD.-^ {\T](k,„ks)f)) -0.000133/: 



2kbvl 

Forz= 1.5,7 = 0.7, Az 



0.2 we then have 



(ir„(^)P) 



1. Ox 10^ (h-^Mpcf' ' 



(39) 



(40) 



where the average is over (ka, ks) at fixed k-. Since this is consid- 
erably larger than P„o, and Plss ^ ^rsi^) must be subtracted from 
the measured r(^) for each k. This can be done by fitting the ob- 
served r(^) with the f{ka,ks)/kr form given by ( 1351 ). In princi- 
ple, this subtraction can be done with a precision cty ~ P„„i/ y/N^ 
where A^^ is the number of radial modes over which the determi- 
nation of r„ is made. If this precision can be reached in practice, 
the subtracted power spectrum is not degraded below the spec- 
trum in the absence of radio sources. 

6. Sensitivity to cosmological parameters 

In Sec. |4] we studied the impact of the various telescope con- 
figurations on power spectrum reconstruction. Fig. |5] shows the 
various power spectra, and it allows us to rank visually the con- 
figurations in terms of electronics noise fraction. The differences 
in P„oi will translate into differing precisions in the reconstruc- 
tion of the BAO peak positions and in the estimation of cosmo- 
logical parameters. 



(34) 6.1. BAO peak precision 

In order to estimate the precision with which BAO peak posi- 
tions can be measured, we have generated power spectra that are 
the sum of the expected HI power spectra and the noise spectra 
calculated in the Section 14.11 The peaks in the generated spec- 
tra were then determined by a fitting procedure and the recon- 
structed peak positions compared with the generated peak posi- 



(35) tions. 
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To this end, we used a metho d similar to the one estabHshed 
in iBlake and Glazebrookl (l2003h . We generated power spectra 
for slices of Universe with a half-sky coverage and a redshift 
depth, Az = 0.2 for < z < 1.6. The power spectrum used in the 
simulation is the sum of HI signal term, coiTesponding to WH 
and noise term derived from ( |26] |. The simulated power spectra 
HI is: 



= 1 + 



Pref(k±,k\i) 

Akexp{-(k/T)")sm 



27T- 



\kl 



k\ 



BAO^ 



Jfc2 
'^BAOW ) 



(41) 



where k - Jk\ + A:?, the parameters A, a and t are adjusted to 

the formula presented in lEisenstein and Hul ( Il998l) . Pref(k±, Ay) 
is the envelop curve of the HI power spectrum without baryonic 
oscillations. The parameters ksAOi. and kgAOW ^^ '^he inverses of 
the oscillation periods in k-space. In the simulations, we used the 
following values for these parameters: A - 1.0, r = 0.1 hMpc^\ 
a = 1.4 and kBAO± - ksAOW = 0.Q6Q hMpc'K 

Each simulation is performed for a given set of param- 
eters which are: the system temperature, r^vj, an observation 
time, Tohs, an average redshift and a redshift depth, Az - 0.2. 
Then, each simulated power spectrum is fitted with a two di- 
mensional normalized function P,„r{kj_, k\\)IPref{k^, k\^\) which is 
the sum of the signal power spectrum of flTb and the normal- 
ized noise power spectrum defined by ( |26T l multiplied by a lin- 
ear term, aq + a^k. The upper limit fc„,aj. in k of the fit cor- 
responds to the approximate position of the linear/non-linear 
transition. This limit is established on the basis of the crite- 
rion discussed in Blake and Glazebrook (2003) . In practice, we 
used for the redshifts z = 0.5, 1.1 and 1.5 respectively k^a^ - 
0.145 hMpc-\ 0.19 hMpc-^ and 0.23 /zMpc"'. 

Figure Q shows the result of the fit for one of theses simu- 
lations. Figure |8] histograms the recovered values of ksAOi. and 
kBAOW fo'" 100 simulations. The widths of the two distributions 
give an estimate the statistical errors. 

In addition, in the fitting procedure, both the parameters 
modeling the signal A, t, a and the parameter correcting the 
noise power spectrum (aq, ai) are floated to take into account 
the possible ignorance of the signal shape and the uncertainties 
in the computation of the noise power spectrum. In this way, we 
can correct possible imperfections and the systematic uncertain- 
ties are directly propagated to statistical errors on the relevant 
parameters kBAO± and ksAOW- ^Y subtracting the fitted noise con- 
tribution to each simulation, the baryonic oscillations are clearly 
observed, for instance, on Fig.|9l 

In our comparison of the various telescope configurations, 
we have considered the following cases for Az - 0.2 slices with 
0<z< 1.6. 

- Simulation without electronics noise: the errors of power 
spectrum are directly related to the number of modes in the 
Az = 0.2 slice of 27rsr O. 

- Unpacked cylinder array (configuration a) 

- Packed cylinder array: (configuration b). 

- Dish array: (configuration c). 

Table [T] summarizes the result. The ranking of the three con- 
figurations is the same as that which we can deduce from the 
noise levels in Fig. |5] the best sensitivity is obtained with the 
packed cylinder array, then with the packed dish array and fi- 
nally with the unpacked cylinder array. However, none of the 
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Fig. 7. ID projection of the power spectrum for one simulation. 
The HI power spectrum is divided by an envelop curve P{k)ief 
coiTesponding to the power spectrum without baryonic oscilla- 
tions. The dots represents one simulation for a "packed" array of 
cylinders with a system temperature, Tj^.j - 50K, an observation 
time, Toi,s = 4month, a solid angle of Insr, an average redshift, 
z = 1 .5 and a redshift depth, Az - 0.2. The soUd line is the result 
of the fit to the data. 
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Fig. 8. Distributions of the reconstructed wavelength ksAOi. and 
ksAOW respectively, perpendicular and parallel to the line of sight 
for simulations as in Fig. |7] The fit by a Gaussian of the dis- 
tribution (solid line) gives the width of the distribution which 
represents the statistical eiTor expected on these parameters. 



three configurations is limited by the cosmic variance with a 
four-month observation. As the errors scales with lITots, con- 
figurations a) and c) could be significantly improved if the ob- 
servation time is increased up to one year 

6.2. Expected sensitivity on wq and Wa 

The observations give the HI power spectrum in angle-angle- 
redshift space rather than in real space. The inverse of the peak 
positions in the observed power spectrum therefore gives the an- 
gular and redshift intervals corresponding to the sonic horizon. 
The peaks in the angular spectrum are proportional to dT{z)/s 
and those in the redshift spectrum to dniz)! s- The quantities dj, 
dn and s all depend on the cosmological parameters. Figure [TOl 
gives the angular and redshift intervals as a function of redshift 
for four cosmological models. The error bars on the lines for 
{Q.m,^a) = (0.27,0.73) correspond to the expected errors on 
the peak positions taken from Table [1] for the four-month runs 
with the packed array. We see that with these uncertainties, the 
data would be able to measure w at better than the 10% level. 
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Table 1. Sensitivity on the measurement of kgAOi. and kBAO\\ as a function of the redshift z for various telescope configurations. I'*' 
row: simulations without noise with pure cosmic variance; 2"'' row: simulations with same electronics noise for a telescope in an 
unpacked cylinder array configuration a); 3* row: simulations with same electronics noise for a telescope in a packed cylinder array 
configuration b); 4* row: simulations with same electronics noise for a telescope in a packed dish array configuration c). 
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Fig. 9. ID projection of the power spectrum averaged over 100 
simulations of the packed cylinder array b. The simulations are 
performed for the following conditions: a system temperature, 
Tsys = 50K, an observation time, Tohs - lyear, a solid angle of 
2nsr, an average redshift, z = 1.5 and a redshift depth, Az - 0.2. 
The HI power spectrum is divided by an envelop curve P{k)ref 
corresponding to the power spectrum without baryonic oscilla- 
tions and the background estimated by a fit is subtracted. The 
errors are the RMS of the 100 distributions for each k bin and 
the dots are the mean of the distribution for each k bin. 
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Fig. 10. The two "Hubble diagrams" for BAO experiments. 
The four falling curves give the angular size of the acoustic 
horizon (left scale) and the four rising curves give the red- 
shift interval of the horizon (right scale). The solid lines are for 
{Q.M, i^A, w) = (0.27, 0.73,-1), the dashed for (1, 0, -1) the dot- 
ted for (0.27,0, -1), and the dash-dotted for (0.27,0.73, -0.9), 
The eiTor bars on the solid curve correspond to the four-month 
run (packed aiTay) of Table [T] 



To estimate the sensitivity to parameters describing dark en- 
ergy equation of state , we f ollow the procedure explained in 
iBlake and Glazebrookl(l2003b . We can introduce the equation of 
state of dark energy, w{z) = wq + w^ • z/(l + z) by replacing Qa 
in the definition ofdjiz) and dH(z), (HJ by: 



Qa 



Q" exp 



1 + w(z') 
1 +z' 



dz' 



(42) 



where Q^ is the present-day dark energy fraction with respect to 
the critical density. Using the relative errors on kBAo± and ksAow 
given in Tab. [U we can compute the Fisher matrix for five cos- 
mological parameter: (Q„,, Q./,, h, wq, Wa)- Then, the combination 
of this BAO Fisher matrix with the Fisher matrix obtained for 
Planck mission, allows us to compute the errors on dark energy 
parameters. The Planck Fisher matrix is obtained for the 8 pa- 
rameters (assuming a flat universe): Q.,„, Q.b, h, wq, Wa, erg, «, 
(spectral index of the primordial power spectrum) and t (opti- 
cal depth to the last-scatter surface). The expected errors and the 
Figure of Merit, the inverse of the area in the 95% confidence 



level contours (see Fig.fTTTi are summarized in Tab.|2] The rank- 
ing between the various configurations in terms of performances 
is consistent with the level of electronics noise observed on the 
power spectra of Fig. |5] 

For an optimized project over a redshift range, 0.4 < z < 1 .6, 
with a total observation time of 4 years, the packed cylinders 
have a precision of 6% on wq and 25% on Wa- Finally, Fig. [12] 
shows a comparison of different BAO projects, with a set of pri- 
ors on (Q„,, Q/,, h) corresponding to the expected precision on 
these parameters in early 2010's. This BAO project based on HI 
intensity mapping is clearly compet itive with the nex t generation 
of optical surveys s uch as SDSS-III (ISDSS-IIIl2008l) or WFMOS 
(iBasset et al.ll2005h . 

7. Conclusion 

In this paper we have discussed the measurement of the HI 
power spectrum with interferometric surveys. By presenting a 
simplified procedure for reconstructing the Fourier space map of 
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Fig. 11. Icr and Icr confidence level contours in the parameter 
plane (wo, Wa) for three configurations : a) unpacked array of 
north-south oriented cylinders (black line), b) packed aiTay of 
north-south oriented cylinders (blue line) and c) packed array of 
dishes (red line). 
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Fig. 12. Icr and 2cr confidence level contours in the parameter 
plane (wo, Wa) for three BAO projects: WFMOS project (blue 
dotted line), SDSS-III (red dashed line) and packed cylinder ar- 
ray with HI intensity mapping (black solid line). 



the GHz sky, we derived expressions for the noise power spec- 
trum. Adding this spectrum to the expected HI power spectrum, 
we determined with what precision the positions of the BAO 
peaks can be measured. This led to a sensitivity to cosmologi- 
cal parameters that is competitive with other BAO projects. 

In calculating the noise power spectrum, we considered three 
interferometer arrays using cylinders or dishes. Packed-cylinder 
and packed-dish arrays have comparable noises so the choice for 
one or the other must be based on cost and technical considera- 
tions, with cylinder arrays having simpler mechanical design but 
requiring more receivers. 
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Table 2. Sensitivity on wo and Wa and Figure of Merit (FoM) for 
various telescope configurations. P' row: simulations without 
noise with pure cosmic variance; 2"'' row: simulations with same 
electronics noise for a telescope in an unpacked cylinder ar- 
ray configuration a); 3* row: simulations with same electronics 
noise for a telescope in a packed cylinder array configuration b); 
4"* row: simulations with same electronics noise for a telescope 
in a packed dish array configuration c); 6* row: Configuration 
b) for an optimized 4-year survey with 0.4 < z < 1.6. 





cr(wo) 


cr(Wa) 


FoM (95%) 


No Noise 

a) Unpacked cylinder array 

b) Packed cylinder array 

c) Packed Dish array 
Packed cylinder array - 4 years 


0.059 
0.17 
0.091 
0.15 
0.063 


0.24 
0.69 
0.36 
0.64 
0.25 


136 
10.1 
69.1 
13.6 
119 



